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ABSTRACT 

This paper presents a sensitivity analysis of the Orthoglide, 
a 3-DOF translational Parallel Kinematic Machine. Two com- 
plementary methods are developed to analyze its sensitivity to its 
dimensional and angular variations. First, a linkage kinematic 
analysis method is used to have a rough idea of the influence 
of the dimensional variations on the location of the end-effector. 
Besides, this method shows that variations in the design parame- 
ters of the same type from one leg to the other have the same in- 
fluence on the end-effector. However, this method does not take 
into account the variations in the parallelograms. Thus, a differ- 
ential vector method is used to study the influence of the dimen- 
sional and angular variations in the parts of the manipulator on 
the position and orientation of the end-effector, and particularly 
the influence of the variations in the parallelograms. 

It turns out that the kinematic isotropic configuration of the 
manipulator is the least sensitive one to its dimensional and an- 
gular variations. On the contrary, the closest configurations to its 
kinematic singular configurations are the most sensitive ones to 
geometrical variations. 

Keywords: Parallel Kinematic Machine, Sensitivity Analysis, 
Kinematic Analysis, Kinematic Singularity. 

NOMENCLATURE 

'Mjj(0,x,y,z) '■ reference coordinate frame centered at O, the in- 
tersection between the directions of the three actuated pris- 
matic joints. 

%j>(P,X ,Y,Z) : coordinate frame attached to the end-effector. 



%i(Ai,Xi,yi,Zi) ■ 
joint. 

P= [px Py Pz\ 

8p 



coordinate frame attached to the i th prismatic 



vector of the Cartesian coordinates of the 



end-effector, expressed in 5^,. 

= [?>p x ?>Py §Pz] r : position error of the end-effector, ex- 
pressed in %},. 

88 = [ 88, 88, 88 z 1 T : orientation error of the end-effector, ex- 



Pi 
8p; 

Li 
a, : 
n : 
di : 
8L, 



pressed in %},. 

displacement of the i th prismatic joint. 

: displacement error of the i' h prismatic joint. 

theoretical length of the i th parallelogram. 

distance between points O and A,-. 

distance between points P and Q. 

theoretical width of the i th parallelogram. 

: variation in the length of the i ,h parallelogram. 



bLij : variation in the length of link ByC,- ; 
dbi : variation in the length of link B/iBq. 
8c, 
8/i : 
8m; 
W; : 
8w, 
8e, 



variation in the length of link Q1Q2. 
parallelism error of links B[\Ba and QiCq. 



: parallelism error of links Bnd\ andB/2C,2- 
direction of links BnCn and Ba.Cn- 



: variation in the direction of links BaCa and BaPa- 
sum of the position errors of points A,-, B,, C,. 
88ai = [88a« 88a;} 88ak] : angular variation in the direction 

of the i th prismatic joint 
88b, = [88b,., 88b„ 88b,,-] : angular variation between B;iB,2 

and the direction of the i th prismatic joint. 
§9o = [ 88ck 88c; y SQciz ] T '■ angular variation between the 



end-effector and Ci\Ca- 



*IRCCyN: UMR n° 6597 CNRS, Ecole Centrale de Nantes, Universite de 
Nantes, Ecole des Mines de Nantes 



8y, = [ 8y, x 8ji y 8y, z ] sum of the orientation errors of the i 



th 
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parallelogram with respect to the / prismatic joint and the 

end-effector. 
DOF : degree-of-freedom. 
PKM : parallel kinematic machine. 

1 Introduction 

For two decades, parallel manipulators have attracted the 
attention of more and more researchers who consider them as 
valuable alternative design for robotic mechanisms. As stated 
by numerous authors, conventional serial kinematic machines 
have already reached their dynamic performance limits, which 
are bounded by high stiffness of the machine components re- 
quired to support sequential joints, links and actuators. Thus, 
while having good operating characteristics (large workspace, 
high flexibility and manoeuvrability), serial manipulators have 
disadvantages of low rigidity and low power. Conversely, par- 
allel kinematic machines (PKM) offer essential advantages over 
their serial counterparts (lower moving masses, higher rigidity 
and payload-to-weight ratio, higher natural frequencies, better 
accuracy, simpler modular mechanical construction, possibility 
to locate actuators on the fixed base). 

However, PKM are not necessarily more accurate than their 
serial counterparts. Indeed, even if the dimensional variations 
can be compensated with PKM, they can also be amplified con- 
trary to with their serial counterparts, [lj. Wang et al. [2] studied 
the effect of manufacturing tolerances on the accuracy of a Stew- 
art platform. Kim et al. used a forward error bound analysis 
to find the error bound of the end-effector of a Stewart platform 
when the error bounds of the joints are given, and an inverse error 
bound analysis to determine those of the joints for the given error 
bound of the end-effector. Kim and Tsai |4] studied the effect of 
misalignment of linear actuators of a 3-DOF translational paral- 
lel manipulator on the motion of its moving platform. Han et al. 
used a kinematic sensitivity analysis method to explain the 
gross motions of a 3-UPU parallel mechanism, and they showed 
that it is highly sensitive to certain minute clearances. Fan et 
al. @ analyzed the sensitivity of the 3-PRS parallel kinematic 
spindle platform of a serial-parallel machine tool 

This paper aims at analyzing the sensitivity of the Orthoglide 
to its dimensional and angular variations. The Orthoglide is a 
three degree-of-freedom (DOF) translational PKM developed by 
Chablat and Wenger Q. A small-scale prototype of this manip- 
ulator was built at IRCCyN. 

Here, the sensitivity of the Orthoglide is studied by means of 
two complementary methods. First, a linkage kinematic analysis 
is used to have a rough idea of the influence of the dimensional 
variations to its end-effector and to show that the variations in 
design parameters of the same type from one leg to another one 
have the same influence on the location of the end-effector. Al- 
though this method is compact, it cannot be used to know the 
influence of the variations in the parallelograms. Thus, a differ- 



ential vector method is used to study the influence of the dimen- 
sional and angular variations in the parts of the manipulator, and 
particularly variations in the parallelograms, on the position and 
the orientation of its end-effector. 

In the isotropic kinematic configuration, the end-effector of 
the manipulator is located at the intersection between the direc- 
tions of its three actuated prismatic joints, and the condition num- 
ber of its kinematic Jacobian matrix is equal to one, [8|. It is 
shown that this configuration is the least sensitive one to geo- 
metrical variations, contrary to the closest configurations to its 
kinematic singular configurations, which are the most sensitive 
to geometrical variations. 

2 Manipulator Geometry 

The kinematic architecture of the Orthoglide is shown in 
FigJTj It consists of three identical parallel chains that are for- 
mally described as PRP a RR, where P, R and P a denote the pris- 
matic, revolute, and parallelogram joints respectively. The mech- 
anism input is made up of three actuated orthogonal prismatic 
joints. The output body (with a tool mounting flange) is con- 
nected to the prismatic joints through a set of three kinematic 
chains. Inside each chain, one parallelogram is used and ori- 
ented in a manner that the output body is restricted to transla- 
tional movements only. 




Figure 1 . Basic kinematic architecture of the Orthoglide 



The small-scale prototype of the Orthoglide is designed such 
that its prescribed performances are Cartesian velocity of 1 .2 m/s 
and an acceleration of 17 m/s 2 . The desired pay load is 4 kg (spin- 
dle, tool, included). The size of its prescribed cubic workspace, 
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C u , is 200 x 200 x 200 mm, where the velocity transmission fac- 
tors are bounded between 1/2 and 2. The three legs are sup- 
posed to be identical. According to [7|, the nominal lengths, Li, 
and widths, di, of the parallelograms, and the nominal distances, 
r,, between points C, and the end-effector P are identical, i.e.: 
L = L\ = L2 = L3 = 310.58 mm, d = d\ = d% = d$ = 80 mm, 
r = r\ = r% = r$ = 3 1 mm. 

As depicted in Figf2] Q\ and Q2, vertices of C„, are de- 
fined at the intersection between the Cartesian workspace bound- 
ary and the axis x = y = z expressed in the reference coordinate 
frame Kj,. Q\ and Q2 are the closest points to the singularity sur- 
faces. Their Cartesian coordinates, expressed in are equal to 
(-73.21,-73.21,-73.21) and (126.79,126.79,126.79), respectively. 




Cartesian 
workspace 



Figure 2. Cartesian workspace, C„, points Q\ and Q2 



The parts of the manipulator are supposed to be rigid-bodies 
and there is no joint clearances. Its i th leg is depicted in Figf3] 
and is made up of one prismatic joint, one parallelogram, and 
three revolute joints, which generate one DOF each. According 
to Karouia et al. Q, the manipulator is isostatic because its legs 
are identical and have five DOF each. Thus, the results obtained 
by the sensitivity analysis methods used in this paper are mean- 
ingful. 



3 Sensitivity Analysis 

Two complementary methods are used to study the sensitiv- 
ity of the Orthoglide. First, a linkage kinematic analysis is used 




Figure 3. Morphology of the i' h leg of the Orthoglide 



to have a rough idea of the influence of the dimensional vari- 
ations to its end-effector. Although this method is compact, it 
cannot be used to know the influence of the variations in the par- 
allelograms. Thus, a differential vector method is used to study 
the influence of the dimensional and angular variations in the 
parts of the manipulator, and particularly variations in the paral- 
lelograms, on the position and the orientation of its end-effector. 

3.1 Linkage Kinematic Analysis 

This method aims at computing the sensitivity coefficients 
of the position of the end-effector, P, to the design parameters 
of the manipulator. First, three implicit functions depicting the 
kinematic of the manipulator are obtained. A relation between 
the variations in the position of P and the variations in the design 
parameters follows from these functions. Finally a sensitivity 
matrix, which gathers the sensitivity coefficients of P, follows 
from the previous relation written in matrix form. 

3.1.1 Formulation Figure Q]depicts the design param- 
eters taken into account. Points A\, A%, and A3 are the bases of 
the prismatic joints. Their Cartesian coordinates, expressed in 
%b, are ai, a 2 , and 83, respectively. 

ai = [-fli 0] r ; a 2 = [0 -a 2 0] T ; a 3 = [00 -a 3 ] T 

where a, is the distance between points A, and O, the origin of 
(Rj,. Points B\, B2 and B3 are the links between the prismatic and 
parallelogram joints. Their Cartesian coordinates, expressed in 
Hb are: 





-a\ +pi 




t>2x 




t>3x 


bi = 


bly 


M = 


-CI2 + P2 


,b 3 = 


Hy 








b 2z 




-fl 3 + P3 



where p, is the displacement of the / prismatic joint. b\ y and 
b\ z are the position errors of point B\ according to y and z axes. 
b2x and b2 z are the position errors of point B2 according to x and 
z axes. b^ x and bj y are the position errors of point B3 according 
to x and y axes. These errors result from the orientation errors 
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of the directions of the prismatic actuated joints. The Cartesian 
coordinates of C\, C 2 , and C3, expressed in are the following: 

ci= [px-rx 0] r ; c 2 = [0 Py -r 2 0] T ; c 3 = [0 p z - r 3 ] T 



where p = [ p x p y p z ] T is the vector of the Cartesian coordinates 
of the end-effector P, expressed in Kj,. 

The expressions of the nominal lengths of the parallelo- 
grams follow from eq.(Q}. 



with 



L 1 -=||c i -b i || 2 , i= 1,2,3 



(1) 



where L, is the nominal length of the / parallelogram and 1 1 - 1 1 2 is 
the Euclidean norm. Three implicit functions follow from eq.fl} 
and are given by the following equations: 



-n + p x + a 1 - p 1 f + (py - biy) 2 + ( Pz — b\ z ) 2 — L 2 =0 



F 2 = ( P x ~ b 2x ) 2 + {-r 2 + Py + a 2 - p 2 ) 2 + (p z - b 2z ) 







^3 = (Px - b ix f + (p y - b 3y ) 2 + (-r 3 +p z + a 3 -p 3 ) 2 -Lj = 

By differentiating functions F\, F 2 , and F 3 , with respect to 
the design parameters of the manipulator and the position of the 
end-effector, we obtain a relation between the positioning error 
of the end-effector, 8p, and the variations in the design parame- 
ters, Sq,. 



bFt = A,8p + B,8q, = , i = 1 ,2,3 



(2) 



with 



A, = [dFi/dp x dFi/dpy dFi/dp z ] 

B, = [dFi/dat dFi/dbiy dFi/db iz dFi/dpt dFi/dL, dF,/d n ] 
8p = [ bp x bp y £>p z ] T 

8q/ = [So,- 8hj 8kj Spj 8I4 8r,] r 

where 8a,, 8/1,, 8kj, 8p,, 5L,-, and 8r,, depict the variations in a,, 
hi, k{, pi, Li, and r,-, respectively with hi = b\ y , k\ = b\ z , h 2 = b 2x , 
k 2 = b 2z , h 3 = bi x , h = biy. 

Integrating the three loops of eq.© together and separating 
the position parameters and design parameters to different sides 
yields the following simplified matrix form: 



A8p + B8q = 



[AfA^fe 



B 



Bi 
B 2 
B 3 



p3x3 



c 3xl8 



Sq = [Sqf Sqj 8q 



)18xl 



Equation (O takes into account the coupling effect of the 
three independent structure loops. According to Q, A is the par- 
allel Jacobian kinematic matrix of the Orthoglide, which does not 
meet parallel kinematic singularities when its end-effector cov- 



ers C„. Therefore, A is not singular and its inverse, A , exists. 
Thus, the positioning error of the end-effector can be computed 
using eq.©. 



where 



A 'B: 



8p = C 8q 



dp x /dai dp x /dh\ 
dp y /dai dpy/dh\ 
dp z /da\ dp z /dh\ 



(4) 



dp x /dr3 
dpy/dri 
dp z /dr$ 



el 3xl8 (5) 



(3) 



represents the sensitivity matrix of the manipulator. The 
terms of C are the sensitivity coefficients of the Cartesian coor- 
dinates of the end-effector to the design parameters and are used 
to analyze the sensitivity of the Orthoglide. 

3.1.2 Results of the Linkage Kinematic Analysis 

The sensitivity matrix C of the manipulator depends on the posi- 
tion of its end-effector. 

Figures |H|5]|6] and [7] depict the mean of the sensitivity co- 
efficients of p x , py, p z , and p, when the end-effector covers C„. 
It appears that the position of the end-effector is very sensitive 
to variations in the position of points A, , variations in the lengths 
of the parallelograms, L;, variations in the lengths of prismatic 
joints, p,, and variations in the position of points C, defined by 
Ti (see Figf5J. However, it is little sensitive to the orientation 
errors of the direction of the prismatic joints, defined by param- 
eters biy,b\ z ,b 2x ,b 2z ,b3, x ,bj, y . Besides, it is noteworthy that p x 
(p y , p z , respectively) is very sensitive to the design parameters 
which make up the 1* (2 nd , 3 rd , respectively) leg of the manipu- 
lator, contrary to the others. That is due to the symmetry of the 
architecture of the manipulator. Henceforth, only the variations 
in the design parameters of the first leg of the manipulator will 
be taken into account. Indeed, the sensitivity of the position of 
the end-effector to the variations in the design parameters of the 
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<dpx/dqc> [|j.m/(j,m] 



U 



L.5T 



0.5 



<dpy/dq>> [fam/fam] 
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Figure 4. Mean of sensitivity of Figure 5. Mean of sensitivity of 
p x throughout C u 
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p-y throughout C ( , 

<cp/dqi> [\xwj\im] 
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Figure 6. Mean of sensitivity of Figure 7. Mean of sensitivity of p 



p z throughout C„ 



throughout C M 



Figure 8. Sensitivity of p x to the Figure 9. Sensitivity of p y to the 



variations in the I s leg 

2 




Ql -50 



50 100 Qi 150 



variations in the 1* leg 




Ql -50 



100 Q, 150 



Figure 10. Sensitivity of p to the Figure 11. Global sensitivity of 



variations in the 1 leg 



p. Px> Py, and p z 



second and the third legs of the manipulator can be deduced from 
the sensitivity of the position of the end-effector to variations in 
the design parameters of the first leg. 

Chablat et al. [7] showed that if the prescribed bounds of 
the velocity transmission factors (the kinematic criteria used to 
dimension the manipulator) are satisfied at Qi and Q2, then these 
bounds are satisfied throughout the prescribed cubic Cartesian 
workspace C„. Q\ and Q2 are then the most critical points of C u , 
whereas O is the most interesting point because it corresponds to 
the isotropic kinematic configuration of the manipulator. Here, 
we assume that if the prescribed bounds of the sensitivity coeffi- 
cients are satisfied at Q\ and Q2, then these bounds are satisfied 
throughout C u . 

Figures [8] and [9] depict the sensitivity coefficients of p x 
and p y to the dimensional variations in the 1" leg, i.e.: 
ai,bi y ,b\ z ,pi,L\,ri, along QiQ2- It appears that these coeffi- 
cients are a minimum in the isotropic configuration, i.e.: P = 0, 
and a maximum when P = Q2, i.e.: the closest configuration to 
the singular one. Figure [10] depicts the sensitivity coefficients of 
p along the diagonal QiQ2- It is noteworthy that all the sensi- 
tivity coefficients are a minimum when P = and a maximum 
when P = Q2- Finally, figure [PT1 depicts the global sensitivities 
of p, p x , p y , and p z to the dimensional variations. It appears that 
they are a minimum when P = 0, and a maximum when P = Qz. 

Figures[T2]and[T3]depict the sensitivity coefficients of p x and 
p in the isotropic configuration. In this configuration, the posi- 
tion error of the end-effector does not depend on the orientation 
errors of the directions of the prismatic joints because the sensi- 
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Figure 12. Sensitivity of p x in Figure 13. Sensitivity of p in the 
the isotropic configuration isotropic configuration 



tivity of the position of P to variations in b\y,b\ z ,b2x,b2zi^x-,b^y 
is null in this configuration. Besides, variations in p x , p y , and 
p z are decoupled in this configuration. Indeed, variarions in p x , 
(p y , p z , respectively) are only due to dimensional variations in 
the l ir , (2 nd , 3 rd , respectively) leg of the manipulator. The corre- 
sponding sensitivity coefficients are equal to 1 . It means that the 
dimensional variations are neither amplified nor compensated in 
the isotropic configuration. 

Figures[l4]and[T5]depict the sensitivity coefficients of p x and 
p when the end-effector hits Q2 (P = Q2)- In this case, variations 
in p x , p y , and p z are coupled. For example, variations in p x are 
due to both dimensional variations in the 1 " leg and variations in 
the 2 nd and the 3 rd legs. Besides, the amplification of the dimen- 
sional variations is important. Indeed, the sensitivity coefficients 
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Figure 16. Variations in O — A; chain 



of p are close to 2 in this configuration. For example, as the sen- 
sitivity coefficient relating to L\ is equal to 1.9, the position error 
of the end-effector will be equal to 19/jm if SLi = 10/um. More- 
over, we noticed numerically that £>2 configuration is the most 
sensitive configuration to dimensional variations of the manipu- 
lator. 

According to figures [4] -[71 EH -Q2] variations in design pa- 
rameters of the same type from one leg the other have the same 
influence on the location of the end-effector. 

However, this linkage kinematic method does not take into 
account variations in the parallelograms, except the variations in 
their global length. Thus, a differential vector method is devel- 
oped below. 



3.2.1 Formulation The schematic drawing of the z 
leg of the Orthoglide depicted in Fig|3] is split in order to de- 
pict the variations in design parameters by a vectoriel form. The 
closed-loop kinematic chains O — A; — B,- — By — Cy — C, —P,i = 
1,2,3, j = 1,2, are depicted by Figs fT6lfT9l "Kj is the coordinate 
frame attached to the i th prismatic joint. o,aj,bj,by,Cy,C;,p, 
are the Cartesian coordinates of points 0,Aj,Bj,Bjj,Cij,Ci,P, re- 
spectively, expressed in Kj and depicted in Figf3] 

According to Fig [16] 

a,--o = R,-(ao + 8a/) (6) 



3.2 Differential Vector Method 

In this section, we perfect a sensitivity analysis method 
of the Orthoglide, which complements the previous one. This 
method is used to analyze the sensitivity of the position and the 
orientation of the end-effector to dimensional and angular vari- 
ations, and particularly to the variations in the parallelograms. 
Moreover, it allows us to distinguish the variations which are re- 
sponsible for the position errors of the end-effector from the ones 
which are responsible for its orientation errors. To develop this 
method, we were inspired by a Huang & al. work on a parallel 
kinematic machine, which is made up of parallelogram joints too 

una. 

First, we express the dimensional and angular variations in 
vectorial form. Then, a relation between the position and the 
orientation errors of the end-effector is obtained from the closed- 
loop kinematic equations. The expressions of the orientation and 
the position errors of the end-effector, with respect to the vari- 
ations in the design parameters, are deduced from this relation. 
Finally, we introduce two sensitivity indices to assess the sen- 
sitivity of the position and the orientation of the end-effector to 
dimensional and angular variations, and particularly to the paral- 
lelism errors of the bars of the parallelograms. 



where arj is the nominal position vector of A, with respect to O 
expressed in 8a, is the positioning error of A,. R, is the trans- 
formation matrix from 'Mj to 5^,. I3 is the (3 x 3) identity matrix 
and 



Ri = h 



R? 



0-1 

1 
0-10 



R, 



1 

00 1 

1 00 



According to FigfTTl 

hi - a,- = R,-(pi 4- 8p,)ei + R/SOa, x (p, + 8p,-)ei (7) 

where p, is the displacement of the i th prismatic joint, Sp,- is its 
displacement error, 80a/ = [SGa/* 58a (v S8a«] is the angular 
variation of its direction, and 





"l" 




"0" 


ei = 





; e 2 = 












1 



m = 



1 if j = 1 

-1 if/ = 2 
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Figure 1 7. Variations in A, — B; chain 




Figure 1 8. Variations in fi, — Bij — Cij chain 



According to FigJTSl 

b, = R,[I 3 + m Ai x]^(j)(d/2 + 8b i /2)[h + 88 B/ x]e$) 
cij -bij = LjWj + SLj/Wj + L/8wj (9) 

where a! is the nominal width of the parallelogram, 8£>, is the 
variation in the length of link B,iZ?;2 and is supposed to be equally 
shared by each side of B,. 88b, = [ 80b,.y 80b, v 80b,, ] T is the 
orientation error of link BuBq with respect to the direction of 
the i th prismatic joint, L,- is the length of the i th parallelogram, 
8L{j is the variation in the length of link ByCy, of which w,- is 
the direction, and Sw,- is the variation in this direction, orthogonal 
to w,-. 

According to Fig [19] 

dj -c i = R i [h + 8Bx]fe(j)(d/2 + 8c i /2)[I 3 + m C iX}e$lQ) 
d p = [I 3 + 86 x ]R, (c + 8ci) (11) 



where 8c, is the variation in the length of link CnCa, which 
is supposed to be equally shared by each side of C,\ 80c; = 




Figure 1 9. Variations in Cy — C, — P chain 



[80c/.v 80c,y 88ci z ] is the orientation error of link C/1Q2 with 
respect to link C,P. Co is the nominal position vector of C\ with 
respect to end-effector P, expressed in %j, 8c,- is the position error 
of Q expressed in %j, and 80 = [ 80 A 80 v 80 ; ] T is the orientation 
error of the end-effector, expressed in 

Implementing linearization of eqs.(lo*l[T0ii and removing the 
components associated with the nominal constrained equation 
po = R;(ao + p,ei - c ) +L,w,-, yields 

8p = p - p 

= R,-(8e, + p,(80 Al - x ei) + 5(y') d l 2 ( §e Ai x e 2 ) + (12) 
%(J) d/2 (SYiXe 2 )+^(j) 8m,/2e 2 ) + 
8L, 7 w,- + L,-8w/ - 80 x R,-(co + d/2 ^(j) e 2 ) 

where 

8p is the position error of the end-effector of the manipulator. 
8e, = 8a, + 8p,ei — 8c,- is the sum of the position errors of points 

A,-, Bi, and C, expressed in Kj. 
Sy,- = 80b, — 80c,- is the sum of the orientation errors of the i' h 

parallelogram with respect to the i th prismatic joint and the 

end-effector. 

8m, = 8/7, — 8c, corresponds to the parallelism error of links 
BaQi and BqPq, which is depicted by Fig|20l 

Equation ( fT2l shows the coupling of the position and orien- 
tation errors of the end-effector. Contrary to the orientation er- 
rors, the position errors can be compensated because the manip- 
ulator is a translational 3-DOF PKM. Thus, it is more important 
to minimize the geometrical variations, which are responsible for 
the orientation errors of the end-effector than the ones, which are 
responsible for its position errors. 

The following equation is obtained by multiplying both 
sides of eq. (fT2l by wf and utilizing the circularity of hybrid 
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c„ 



"C 

+ 



B„ 



L+bL n 



nominal 
parallelogram 



L,+8L n 



81, 



\ 8m, 



real parallelogram 



•th 

Figure 20. Variations in the v parallelogram 



product. 



The determinant of D will be null if the normal vectors to the 
plans, which contain the three parallelograms respectively, are 
collinear, or if one parallelogram is flat. Here, this determinant 
is not null when P covers C„ because of the geometry of the ma- 
nipulator. Therefore, D is nonsingular and its inverse D 1 exists. 

As (R,e2 x w,') r _LR,e2, 50A iZ and 5y iZ , the third components 
of S8ai and 8y expressed in Kj, have no effect on the orientation 
of the end-effector. Thus, matrix Jee can be simplified by 
eliminating its columns associated with 80a ( v and Sy,v, i = 1,2,3. 
Finally, eighteen variations: 8/,-, 8m,, 80a«, S0a!>, 8y„, 8y n , 
i = 1,2,3, should be responsible for the orientation error of the 
end-effector. 

Position Error Mapping Function: By addition of eqs.(fT3l 

written for j = 1 and j = 2, and for the i th kinematic chain, a 
relation is obtained between the position error of the end-effector 
and the variations in design parameters, which does not depend 
on 8y . 



wf8p - wfRSc + p^Rid x w i ) r R i 80Ai + ^(;) d/2 (13) 
(R,e 2 x w ( -) r Ri(S0Ai + H)+W) Smi/2 wfR,e 2 
+8L,j - (Ri(co + W) d/2 e 2 ) x w.-fse 



Orientation Error Mapping Function: By substraction of 
eqs.([T3l written for j = 1 and j = 2, and for the i' h kinematic 
chain, a relation is obtained between the orientation error of the 
end-effector and the variations in design parameters, which is 
independent of the position error of the end-effector. 

c/(R,e2 x w,) r 80 = 8/, +c/(R,e2 x w ! ) r R,(80A, + 8y,) + 8m,- wf R,e2 

(14) 



where 8Z, = 8L,i — 8L,2, the relative length error of links B/iCn 



and BiiCu, depicts the parallelism error of links BuBq and Qi C,2 
as shown in Figf20] Equation ( fT~4-b can be written in matrix form: 



80 — Jee^e 



(15) 



with 



Jee = D E 



D = d 





"Ei " 


; E = 


. . . E 2 • • • 




E 3 



"(Rie 2 xwi) r 
(R 2 e 2 x w 2 ) r 
_(R 3 e 2 x W3) r _ 

E,- = [1 wfR ie2 d(R,e 2 x w,) r R, d(R,e 2 x w,) r R,] 

80 is the orientation error of the end-effector expressed in Kj,, 
and £e = (z^.z^.z^) 7 such that ee; = (8Z,,8m;,80j,,8yf ) T . 



wf8p = 8L; + wf R;8e ; + p,(R,ei x Wi) T R/80A; - (R;c x w,) r 80 

(16) 

Equation ( [Tol l can be written in matrix form: 



8p = 3, 



with 



F 
G 
H 

Z n 



F G 



■ JpeEe — [JppJpe] 



Jpe = F 'HJee 



£e 



(17) 



[W1W2W3J 

diagCG;) ; Gi = [ 1 wf R; p^Rpi x w,) r R, ] 
- [RlCo X Wl R2C0 x w 2 R3C0 x w 3 ] 

{z T lA ,z T p2 ,z T pi )\ 8 pi -(8L„8ef,80l,.) r 



8L, = (8L,i +8Lq)/2 is the mean value of the variations in 
links BuCii and BaCa, i.e.: the variation in the length of the 
f h parallelogram. z p is the set of the variations in design pa- 
rameters, which should be responsible for the position errors of 
the end-effector, except the ones which should be responsible for 
its orientation errors, i.e.: Zq. Z p is made up of three kinds of 
errors: the variation in the length of the i th parallelogram, i.e.: 
SL,-,! = 1,2,3, the position errors of points A,, B,, and Q, i.e.: 
Se,-,i = 1,2,3, and the orientation errors of the directions of the 
prismatic joints, i.e.: 80^/, i= 1,2,3. Besides, F is nonsingular 
and its inverse F~' exists because F corresponds to the Jacobian 
kinematic matrix of the manipulator, which is not singular when 
P covers C u , 0- 

According to eq.(fT6ll and as (R,ei x w/) _LR;ei, matrix J pp 
can be simplified by eliminating its columns associated with 
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SQaut, i = 1,2,3. Finally, thirty-three variations: 8L,-, 8e, x , Se,-,,, 
Se,>., 80Ai.r, §8a!>, 80a,v, 8/,-, 8m,-, Sy,*, Sy^, i = 1,2,3, should be 
responsible for the position error of the end-effector. 

Rearranging matices J pp and J p e, the position error of the 
end-effector can be expressed as: 



8p = Je q = [Ji h h] (£, 



ql E q2 £ c/3 



T \T 



(18) 



with Zqi = (8L,, Se,*, 5e iy , Se,v, S6a«, 86^,, 
andJeR 3 * 33 . 



,80A !z ,8/,,8m,,8Y«,8Yh 



Sensitivity Indices: In order to investigate the influence of 
the design parameters errors on the position and the orientation 
of the end-effector, sensitivity indices are required. According 
to section [3. 1.21 variations in the design parameters of the same 
type from one leg to the other have the same influence on the 
location of the end-effector. Thus, assuming that variations in the 
design parameters are independent, the sensitivity of the position 
of the end-effector to the variations in the design parameter 
responsible for its position error, i.e.: E q n t 2,3)k' is called an d is 
defined by eq.([T9b. 



EE', 



2 k 

imk ' ' 



,11 



(19) 



Likewise, the sensitivity of the orientation of the end- 
effector to the variations in the r 11 ' design parameter responsible 
for its orientation error, i.e.: £9(1,2,3),-, is called v r and follows 
from eq.dT3T>. 



v r = 



\ 



EE4 

;=0<=1 



•i{6j+r) 



r=l,-,6 



(20) 



Finally, ^ can be employed as a sensitivity index of the po- 
sition of the end-effector to the k th design parameter responsible 
for the position error. Likewise, v r can be employed as a sensi- 
tivity index of the orientation of the end-effector to the r** design 
parameter responsible for the orientation error. It is noteworthy 
that these sensitivity indices depend on the location of the end- 
effector. 

3.2.2 Results of the Differential Vector Method 

The sensitivity indices defined by eqs.([T9l) and < [2Qb are used to 
evaluate the sensitivity of the position and orientation of the end- 
effector to variations in design parameters, particularly to varia- 
tions in the parallelograms. 




Figure 21. Sensitivity of the position of the end-effector along Q1Q2, 
(a): to dimensional variations, (b): to angular variations 



Figures l2TT a-b) depict the sensitivity of the position of the 
end-effector along the diagonal Q1Q2 of C u , to dimensional 
variations and angular variations, respectively. According to 
FigBTTa). the position of the end-effector is very sensitive to 
variations in the lengths of the parallelograms, SL,-, and to the 
position errors of points A,, B,, and C, along axis x, of Kj, i.e: 
8e/ A . Conversely, the influence of 8/, and 8m,, the parallelism 
errors of the parallelograms, is low and even negligible in the 
kinematic isotropic configuration. According to FigfJTJb), the 
orientation errors of the prismatic joints depicted by S8a, v and 
89a,,- are the most influential angular errors on the position of 
the end-effector. Besides, the position of the end-effector is not 
sensitive to angular variations in the isotropic configuration. 




Figure 22. Sensitivity of the orientation of the end-effector along Q1Q2, 
(a): to dimensional variations, (b): to angular variations 



Figures l22l a-b) depict the sensitivity of the orientation of 
the end effector, along Q1Q2, to dimensional and angular varia- 
tions. According to Fig l22l a). 8/, and 8m, are the only dimen- 
sional variations, which are responsible for the orientation error 
of the end-effector. However, the influence of the parallelism er- 
ror of the small sides of the parallelograms, depicted by 8/,-, is 
more important than the one of the parallelism error of their long 
sides, depicted by 8m, . 
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According to figures [21] and [22] the sensitivity of the posi- 
tion and the orientation of the end-effector is generally null in 
the kinematic isotropic configuration (p = 0), and is a maximum 
when the manipulator is close to a kinematic singular configura- 
tion, i.e.: P = Q2- Indeed, only two kinds of design parameters 
variations are responsible for the variations in the position of the 
end-effector in the isotropic configuration: 8L, and 8e u . Like- 
wise, two kinds of variations are responsible for the variations in 
its orientation in this configuration: 8/,-, the parallelism error of 
the small sides of the parallelograms, 884™ and 8y, v . Moreover, 
the sensitivities of the pose (position and orientation) of the end- 
effector to these variations are a minimum in this configuration, 
except for 8/,. On the contrary, Q2 configuration, i.e.: P = Q2, is 
the most sensitive configuration of the manipulator to variations 
in its design parameters. Indeed, as depicted by Figsl2Tlandl22l 
variations in the pose of the end-effector depend on all the design 
parameters variations and are a maximum in this configuration. 

Moreover, figures [21] and [22] can be used to compute the 
variations in the position and orientation of the end-effector with 
knowledge of the amount of variations in the design parameters. 
For instance, let us assume that the parallelism error of the small 
sides of the parallelograms, 8/,, is equal to 10/jm. According 
to Figf22]a), the position error of the end-effector will be equal 
about to 3/jm in Q\ configuration (P = Q\). Likewise, according 
to FigfJTJb), if the orientation error of the direction of the f h pris- 
matic joint round axis yi of %i is equal to 1 degree, i.e.: 80^, v = 1 
degree, the position error of the end-effector will be equal about 
to 4.8 mm in Q2 configuration. 

4 Conclusions 

In this paper, two methods were developed to analyze the 
sensitivity of the Orthoglide. The first method was used to have a 
rough idea of the influence of variations in the design parameters 
on the location of the end-effector. It showed that variations in 
design parameters of the same type from one leg to the other have 
the same influence on the end-effector. The second method was 
used to study the sensitivity of the position and the orientation of 
the end-effector to the dimensional and angular variations, and 
particularly to variations in the parallelograms. 

It turns out that the kinematic isotropic configuration of the 
manipulator is the least sensitive one to its dimensional and an- 
gular variations, contrary to the closest configurations to its kine- 
matic singular configurations. Besides, the position of the end- 
effector is very sensitive to variations in the lengths of the par- 
allelograms and variations in the lengths of the prismatic joints. 
On the contrary, the parallelism errors of the bars of parallelo- 
grams are little influential on the position of the end-effector. As 
far as the orientation of the end-effector is concerned, it is more 
sensitive to the parallelism errors of the small bars of the paral- 
lelograms than to the ones of its long bars. 

Therefore, these results should help the designer of the Or- 



thoglide to synthesize its dimensional and angular tolerances. 
Moreover, they can be used to study the robustness of the ma- 
nipulator. Finally, one of the next steps in our research work is to 
study the relation between the sensitivity and the rigidity of the 
manipulator. 
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